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j^ Abstract. Measurements of London penetration depth are a sensitive tool to study 

S multi-band superconductivity and it has provided several important insights to the 

behavior of Fe-based superconductors. We first briefly review the "experimentalist 
- friendly" self-consistent Eilenberger model that relates the measurable superfluid 
Q density and structure of superconducting gaps. Then we focus on the BaFe2As2 - 

O derived materials, for which the results are consistent with 1) two distinct 

superconducting gaps; 2) development of strong in-plane gap anisotropy with the 
^~~^ departure from the optimal doping; 3) appearance of gap nodes along the c— direction 

J^ in a highly overdoped regime; 4) significant pair-breaking, presumably due to charge 

r — doping; 5) fully gapped (exponential) intrinsic behavior at the optimal doping 

\C if scattering is removed (probed in the "self-doped" stoichiometric LiFeAs); 6) 

^^ competition between magnetically ordered state and superconductivity, which do 

t^ coexist in the underdoped compounds. Overall, it appears that while there are common 

trends in the behavior of Fe-based superconductors, the gap structure is non-universal 
and is sensitive to the doping level. It is plausible that the rich variety of possible gap 
structures within the general s± framework is responsible for the observed behavior. 



1. Introduction 

After the initial discoveries of superconductivity in Fe-based superconductors, first 
in LaFePO with T^ ?« 4 K [1] and then in LaFeAs(Oi_:,F:,) with Tc ^ 26 K [2], T^ as 
high as 55 K has been reported in SniFeAs(0i_2;Fa;) [5]. These findings have triggered 
an intense research aimed at understanding the fundamental physics governing this 
new family of superconductors. Since the original discoveries, at least seven different 
classes of iron-based superconductors have been identified, all exhibiting unconventional 
physical properties [H |5l El [7] . Among those, perhaps the most diverse family is based on 
(AE)Fe2As2 parent compounds (where AE=alkaline earth, denoted the "122" system) 
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in which electron |8], |9l UHl [H], hole [111 [13] and isovalent [m [I5] dopants induce a 
superconductivity "dome" with maximum Tc achieved at fairly high concentrations of 
the dopant (up to 45 % in case of BaK122) and showing coexistence of superconductivity 
and long-range magnetic order on the underdoped side. 

Three plus years into the intense study of the Fe-based superconductors, several 
experimental conclusions supported by a large number of reports can be made. Among 
those, relevant to our discussion of London penetration depth, are: 1) two distinct 
superconducting gaps of the magnitude ratio of about 2; 2) power-law behavior 
of thermodynamic quantities at low temperatures presumably due to pair-breaking 
scattering and gap anisotropy; 3) doping dependent three-dimensional gap structure 
possibly with nodes. For some reviews of various experiments and theories the reader is 
referred to [H El El [ISl [HI [IHl [HI [IH] • Here we focus on London penetration depth as 
a sensitive tool to study superconducting gap structure, albeit averaged over the Fermi 
surface |20]. 



2. Measurements of London Penetration Depth 

There are several ways to measure magnetic penetration depth in superconductors 
[20] . These include muon-spin rotation (/iSR) [2T1 [221 ES], frequency dependent 
conductivity [211 I2S], magnetic force and SQUID microscopy [2S1 [2Z|, measurements of 
the first critical field using either global [281 EH] or local probes [301 EI] , microwave cavity 
perturbation [321 |33l jSH [351 ES], mutual inductance (especially suitable for thin films) 
[37] and self-oscillating tunnel-diode resonator (TDR) [2Sl[Mll3Hllini[lIl[121llSl[31]. 
Each technique has its advantages and disadvantages and only by a combination 
of several different measurements performed on different samples and at different 
experimental conditions, one can obtain a more or less objective picture. In the case of 
Fe-based superconductors, most of the reports agree in general, but still differ on the 
details. Most difficulties come from the necessity to resolve very small variations of the 
penetration depth at the lowest temperatures (below 1 K) while maintaining extreme 
temperature resolution and the ability to measure statistically sufficient number of data 
points over the full temperature range. In addition, it is important to look at a general 
picture by measuring many samples of each composition and over extended doping 
regime. This review is based on the results obtained with a self-oscillating tunnel-diode 
resonator (TDR) that has the advantage of providing stable resolution of about 1 part 
per 10^ [20]. Indeed, comparison of TDR with other techniques, especially bulk probes, 
such as thermal conductivity, is important to come up with a consensus regarding the 
gap structure in pnictide superconductors 

2.1. Tunnel - diode resonator (TDR) 
2.1.1. Frequency - domain measurements: 
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When a non-magnetic conducting sample is inserted into a coil of a tank circuit of 
quality factor Q and resonating at a frequency /o = l/2^T^/L^, it causes changes in 
the resonant frequency and the quality factor. For a flat slab of width 2w and volume 
Vs in parallel field [47j : 

A/ _ /o - / (a) _ Vs [^ ^ tanh(aw) \ 



/o /o 2V; V "^ 

where a = (1 — i)/5 in a normal metal and a = 1/A in a superconductor. Here, (5 is the 
normal skin depth and A is the London penetration depth. Vc is volume of the inductor. 
Solving explicitly, we have for a normal metal: 



^ o T. 2w , „^„ 2w I K"^) 



Af _ Vs f 6 sinh ^ + sin ^'" 
~Y~2Vc\ ~ 2^ cosh ^ + cos ^ 

and, in a superconductor, we obtain the known "infinite slab" solution: 

which allows us to measure both superconducting penetration depth and normal - state 
skin depth (which gives the contact-less measure of resistivity). 

In the case of a finite sample with magnetic susceptibility Xy we can still write: 
A/ = —4:TTx(T)G, where G ~ foVs/2Vc (1 — A^) is a geometric calibration constant that 
includes the effective demagnetization factor A^ [48J . Constant G is measured directly by 
pulling the sample out of the coil [IHl ED] . The susceptibility in the Meissner state can 
be written in terms of London penetration depth A(T) and normal state paramagnetic 
permeability (if local - moment magnetic impurities are present), fJ^iT), as 



where R is the characteristic sample dimension. For a most typical in the experiment 
slab geometry with 26 > 2a planar dimensions and thickness 2d (magnetic field 
excitation is along the d— side), R is approximately given by 



w 
R 



^ ,^ f2dY\ /w\ 2d 



(6) 



with w ~ ah /{a + h). 

Unlike microwave cavity perturbation technique that requires scanning the 
frequency |17], TDR is a self-oscillating resonator always locked onto its resonant 
frequency |19]. A sample to be studied is mounted on a sapphire rod and inserted into 
the inductor coil of the tank circuit. Throughout the measurement the temperature of 
the circuit (and of the coil) is stabilized at ±1 mK. This is essential for the stability in 
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the resonant frequency, which is resolved to about 0.01 Hz. This translates to the ability 
to detect changes in A(T) in the range of a few Angstrom. The ac magnetic excitation 
field, Hac, in the coil is about 20 mOe, which is small enough to ensure that no vortices 
are created and London penetration depth is measured. 

2.1.2. Measurements of the absolute value of \(T): 

The described TDR technique provides precise measurements of the variation of the 
penetration depth, AA(T), but not the absolute value for the reasons described in detail 
elsewhere [50]. However, the TDR technique can be extended to obtain the absolute 
value of the penetration depth, A(T). The idea is to coat the entire surface of the 
superconductor under study with a thin film of a conventional superconductor with 
lower Tc and a known value of A(0). In this work we used aluminum, which is convenient, 
since its T/' ^ 1.2 K is quite low for most of the discussed materials, so it is possible to 
extrapolate to T = and obtain A(0). While the Al film is superconducting it screens 
the magnetic field and the effective penetration depth is [39] : 



A(T) + AAz(T)tanh, , . 

Ae//(T) = Xai{T) y^, (7) 

A^K3^) + A(r)tanh(^j 

However, when it becomes normal, Al layer does not cause any screening because its 
thickness, t, is much less than the normal state skin depth at the TDR operating 
frequency of 14 MHz, where 6ai ~ 75 /im for /9^'=10 fiQ-cra [51]. By measuring the 
frequency shift upon warming from the base temperature, T = Tmin, to T = T^^ where 
Xai -> oo and Eq. Q gives XefjiTf) = t + X{Tf), we can calculate X{Tf). We also 
note that if aluminum coating is damaged (i.e., cracks and/or holes), the result will 
underestimate A(0). This method, therefore, provides the lower boundary estimate. 

2.1.3. Out-of-plane Penetration Depth: 

Finally, the TDR technique can be used to measure out-of-plane component of the 
penetration depth, Ac. For the excitation field Hac \\ c, screening currents flow only in 
the a6-plane and A/ is only related to the in-plane penetration depth, AAa^. However, 
when the magnetic field is applied along the a6-plane {Hac || cib), screening currents flow 
both in the plane and between the planes, along the c-axis. In this case. A/-*- contains 
contributions from both AAat and AAc. For a rectangular sample of thicknesses 2t, 
width 2w and length /, A/-*- is approximately given by Eq. (pi) 



^/ ^ AAa^, AAc _ AAmJi /Q\ 

'Kf^^^ ~^~ Rt ^^ 

where Rb is the effective dimension that takes into account finite size effects [48], Eq. [6] 
Knowing AXab from the measurements with Hac \\ c and the sample dimensions, one can 
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obtain A Ac from Eq. rtSj). However, because 2w > 4 x 2t in most cases, A/ is in general 
dominated by the contribution from AX^b- The subtraction of AAc becomes therefore 
prone to large errors. The alternative, more accurate approach is to measure the sample 
twice |52|. After the first measurement with the field along the longest side / {Hac \\ I), 
the sample is cut along this / direction in two halves, so that the width (originally 2w) 
is reduced to w. Since the thickness 2t remains the same, we can now use Eq. ([s]) to 
calculate AAc without knowing AXab- Note that the length / and width w are in the 
crystallographic a6— plane, whereas the thickness 2t is measured along the c— axis. In 
our experiments, both approaches to estimate AAc(r) produced similar temperature 
dependence, but the former technique had a larger data scatter, as expected. 

3. London Penetration Depth and Superconducting Gap 

In order to describe London penetration depth in multi-gap superconductors and to 
take into account pair-breaking scattering, we use the weak coupling Eilenberger quasi- 
classical formulation of the superconductivity theory that holds for a general anisotropic 
Fermi surface and for any gap symmetry [53]. This method suitable for the analysis of 
the experimental data is described in a few publications of one of us [Ml ESI ESI SD] ; here 
we briefly summarize the scheme. 

In the clean case the Eilenberger equations read: 

vUf = 2Ag/h-2iof, (9) 

-vU*r = 2A*g/h-2iof+, (10) 

(11) 



/ = 


1- 


LOB 


A{r. 


,v) 


= 27rTN{0) ^ 

aj>0 



V{v,v')f{v',r,co))^^, (12) 

-47r|e|A^(0)T Im ^ /vg\ . (13) 



a;>0 

Here, v is the Fermi velocity, LI = V + 27iiA/(f)Q, (pQ is the flux quantum. A(r) is 
the gap function (the order parameter) which may depend on the position kp at the 
Fermi surface (or on v). Eilenberger functions f{r,v,uj), f~^, and g originate from 
Gor'kov Green's functions integrated over the energy variable near the Fermi surface 
to exclude fast spatial oscillations on the scale l/kp; instead f,g vary on the relevant 
for superconductivity scale of the coherence length ^. Functions /, /^ describe the 
superconducting condensate, whereas g represents normal excitations. A^(0) is the total 
density of states at the Fermi level per spin. The Matsubara frequencies are defined by 
hw = 7iT{2n + 1) with an integer n, and ud being the Debye frequency (for phonon- 
mediated superconductivity or a relevant energy scale for other mechanisms). The 
averages over the Fermi surface weighted with the local density of states oc l/\v\ are 
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defined as 

X)= [ ^ /^^, X. (14) 



The order parameter A is related to / in the self-consistency equation Eq. (12). 
Often, the effective couphng V is assumed to be factorizable [57], V{v,v') = 
Vofl{v) fi{v'); this assumption is not always justifiable but it makes the algebra much 
simpler. One looks for the order parameter in the form A{r,T;v) = '${r,T)Q{v). 
Then, the self-consistency Eq. (12) takes the form: 

vl/(r, T) = 27rTiV(0)yo Yl (^(^)/(^' ^' ^)) ■ (15) 

a;>0 

The function Q{v) (or Q{kp)), which describes the variation of A along the Fermi 
surface, is conveniently normalized [58] : 

n^) = i. (16) 



In the absence of currents and fields 11 = 0, and the Eilenberger equations give for 
the uniform ground state: 

/o = /o^ = ^, ^o = y, P' = Al + hW; (17) 

Note that in general, both Aq = '^o(T)Q{kp) and (3 depend on the position kp at the 
Fermi surface. 

Instead of dealing with the effective microscopic electron-electron interaction V, 
one can use in this formal scheme the measurable critical temperature T^ utilizing the 
identity 

-1 rp ^D 



In 



H-2.r^-, (18) 



which is equivalent to the famous relation A(0) = ttTcC"'*' = 2/ia;Dexp(— l/A^(0)Vo) 



7 = 0.577 is the Euler constant. Substitute Eq. ( pL8| ) in Eq. (15) and replace ud with 
infinity due to fast convergence: 



2nT T ^^ \huj 

uj>0 



Now equation for ^(T) reads: 



1 , T, ^ / 1 / fi2 

In 



27rT T ^ V^^ \V^^WTn^' 

a;>0 ^ ^ 
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3.1. London Penetration Depth 

Within microscopic theory, the penetration of weak magnetic fields into 
superconductors is treated perturbatively. Weak supercurrents and fields leave the order 
parameter modulus unchanged, but cause the condensate, i.e., A and the amplitudes 
/ to acquire an overall phase 0{r). Using the method of perturbations, one obtains 
corrections to /cS'o among which we need only 

gi = th^^v-{Ve + 2T,AI<PQ). (21) 



Substituting this in the general expression (13) for the current density, one obtains 
the London relation between the current and the "gauge invariant vector potential" 
a = (j)oVe/2TT + A: 

47rjVc = -(A2)-V. (22) 



Then, the general expression ( 13 ) for the current gives the inverse tensor of the squared 
penetration depth: 

i>^")ik = -^ mo) J2 (^) • (23) 

This result holds at any temperature for clean materials with arbitrary Fermi surface 
and order parameter anisotropies [31]. The temperature dependence of Aq = ^Pfi and 



(3 should be obtained solving Eq. (20). 



Thus, the general scheme for evaluation of A(T) in clean superconductors consists 
of two major steps: first evaluate the order parameter Ao(T) in uniform zero-field state 



for a given gap anisotropy ^{v), then use Eq. (23) with a proper averaging over the 
Fermi surface. The sum over Matsubara frequencies is rapidly convergent and is easily 
done numerically (except in a few limiting situations for which analytic evaluation is 
possible). 

3.1.1. Isotropic A on a general Fermi surface: For majority of materials with 
electron-phonon interaction responsible for superconductivity, the relevant phonons have 
frequencies on the order of uu- One can see that in exchanging such phonons, the 
electron momentum transfer is of the order of hkp', this leads to considerable smoothing 
of A(fcir) [59]. Indeed, in such materials A(A;^) is nearly constant along the Fermi 
surface, the strong possible anisotropy of the latter notwithstanding. In this case, 
commonly called "s— wave" , A is taken as a constant at the Fermi surface: 

(A'),-^ «'<"-"^> 2.rA^^-i. (24) 

LlJ 

We obtain in the Ginzburg-Landau domain: 

{X%\T ^ n) = 2{X%\0) (1-T/T,), t = T/r,. (25) 
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^ = ^ _ pim e-A<o./- « 1.T64 - XmV-te-""" . (26) 
The low-temperature behavior of the penetration depth is given by: 



i^)^ = i^\,'iO) I 1 - 2^/^e-^(°)/^ ) . (27) 

3.1.2. 2D d—wave: As an example of anisotropic A, let us take a relatively simple but 
important case of a (i— wave order parameter on the two-dimensional cylindrical Fermi 
surface: i7 = i7o cos 2ip where if is the properly chosen azimuth angle on the Fermi 
cylinder. With this choice, the gap nodes are at y9 = ±7r/4, ±37r/4. The normalization 



Eq.(16) gives VLq = a/2, so that Aq = \E'v^cos2(y9. The order parameter at T = is now 
given by: 

27r 

2.139. (28) 



A, 


naM 




Tc 


and at T < T^: 




A, 


nax {-'' ) 



g7+0.5 



2.139- 0.927 1^ (29) 

After averaging over Fermi cylinder we obtain for the in-plane penetration depth: 

A-^(T) = A-'^(0)(l-^£-) , (30) 

WiV(OK^ " 

3.2. Eilenberger two-gap Scheme: The 7 - model 

The full-blown microscopic approach based on the Eliashberg theory is quite involved 
and not easy for the data analysis [6OI [61], [621 [63]. Hence, the need for a relatively 
simple but justifiable, self-consistent and effective scheme experimentalists could employ. 
The weak-coupling model is such a scheme. Over the years, the weak-coupling theory 
had proven to describe well multitude of superconducting phenomena. Similar to the 
weak coupling is the "renormalized BCS" model [63] that incorporates the Eliashberg 
corrections in the effective coupling constants. We call our approach a "weak- coupling 
two-band scheme" and refer the reader to original papers where it is clarified that the 
applicability of the model for the analysis of the superfiuid density and specific heat 
data is broader than the traditional weak coupling [M] [5S] . 

The Eilenberger approach can be used to describe self-consistently two-gap 
situation, in which 

A(k) = Ai,2, kGFi,2, (32) 
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where Fi,F2 are two separate sheets of the Fermi surface. Denoting the densities of 
states on the two parts as A'^1,2, we have for a quantity X constant at each Fermi sheet: 

(X) = {X^N^ + X2N2)/N{0) = n,X^ + n^X^ , (33) 

where ni^2 = Ni^2lN{Q)] clearly, ni + n2 = 1. 

The self-consistency equation ( [l2| now takes the form: 

Aj, = 2'nT ^ n^K^f^ = ^ n^A^^A^ ^ — — , (34) 

where z/ = 1,2 is the band index and X^^j. = N{0)V{i', ^) are dimensionless effective 
interaction constants. Note that the notation commonly used in literature for Xl^ 
differs from ours: Xlli = n^Xi,^. 

Turning to the evaluation of Aj,(T), we note that the sum over u in Eq. (34) is 
logarithmically divergent. To deal with this difficulty, we employ Eilenberger's idea of 
replacing fuvD with the measurable Tc. Introducing dimensionless quantities 

^^ ~ 27CT ~ T. 2nt ' ^^^^ 



with t = T/Tc, we obtain: 



1^=1,2 ^^ ^ / 



where A is defined as: 

1.76Te = 2;ttc;Dexp(-l/A), (37) 

or 

1 , T ^2nT , , 

- = In — + > —— . 38 

In terms of X^^, X is expressed as: 

^ ^ 2^1^2(AllA22 - AJ2) /ggN 

niXii + n2X22 - ^/{niXii - n2X22Y - 4nin2Ai2 
For the given coupling constants X^^ and densities of states njy, this system can 
be solved numerically for 5y and therefore provide the gaps Aj, = 2TiT5y{t). Example 
calculations are shown in Fig. [ij First graph in the top row is calculated assuming 
no interband coupling. Naturally, we obtain material with two different transition 
temperatures. The second graph shows the gaps with A12 = 0.05, which features a 
single Tc and quite non-single-BCS-gap temperature dependence of the smaller gap. 
The single BCS gap is shown for comparison. 
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Figure 1. Top row: self-consistent gaps calculated for the indicated X,y^ in a 
compensated metal with /i = 1. First graph is for no inter-band pairing, A12 = 0. 
The second for A12 = 0.05. Bottom row: corresponding partial and total superfluid 
densities. 



Superfluid density: Having formulated the way to evaluate A(T), we turn to the 
London penetration depth given for general anisotropies of the Fermi surface and of A 



by Eq. ( 23 ) ^54j • We consider here only the case of currents in the ab plain of uniaxial or 
cubic materials with two separate Fermi surface sheets, for which the superfluid density, 
P = Al,(0)/A^,(T),is: 



P =7Pi + (l-7)P2, 

00 



p>y 



7 



n=0 



riivf + 77,2^1 



(40) 



where v^ are the averages of the in-plane Fermi velocities over the corresponding band. 
With the discovery of two-gap superconductivity in a number of materials, including 
MgBs [Ml ES], NbaSe [52], VgSi [66], LuaFegSis [67] and ZrBi2 [68] one of the most 
popular approaches to analyze the experimental results has been the so called "a-model" 
[6^ . Developed originally to renormalize a single weak-coupling BCS gap to account 
for strong - coupling corrections [SS], it was used to introduce two gaps, A12, each 
having a BCS temperature dependence, but different amplitudes [Mj- This allowed 
for a simple way to fit the data on the specific heat [64j and the superfluid density, 
p = xpi + (1 — x)p2 [651 EDj. Here, pi,2 are evaluated with Ai^2 = («i,2/1-76)Abcs'(^) 
and X takes into account the relative band contributions. The fitting is usually quite 
good (thanks to a smooth and relatively "featureless" p(T) and the parameters a where 
found to be one larger and one smaller than unity (unless they are both equal to one in 
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the single-gap limit). Although the a-model had played an important and timely role 
in providing convincing evidence for the two-gap superconductivity, it is intrinsically 
inconsistent as applied to the full temperature range. The problem is that one cannot 
assume a priory temperature dependencies for the gaps in the presence of however weak 
interband coupling (required to have single T^). In unlikely situation of zero interband 
coupling, two gaps would have single-gap BCS-like T— dependencies, but will have two 
different transition temperatures. The formal similarity in terms of additive partial 
superfluid densities prompted to name our scheme the "7-model". We note, however, 
that these models are quite different: our 7 that determines partial contributions from 
each band is not just a partial density of states Ui of the a-model, instead it involves 
the band's Fermi velocities. The gaps, Ai 2(T), are calculated self-consistently during 
the fitting procedure. 




Figure 2. Symbols: superfluid density, Ps{T) calculated with A(0) = 200 nm. Solid 
lines represent the fit to a two-gap 7— model, ps — "fpi + (1 — l)p2- Dashed line is 
a single-gap BCS solution. Upper inset: superconducting gaps, Ai(T) and A2(T) 
calculated self-consistently during the fitting. Lower inset: A1/A2 as a function of 
temperature. 



The 7-models can be simplified for a compensated metal, such as clean 
stoichiometric superconductor LiFeAs jH]. To reduce the number of fitting parameters, 
yet capturing compensated multiband structure, we consider a simplest model of two 
cylindrical bands with the mass ratio, /i = mi/m2, whence the partial density of states 
of the first band, ni = /i/(l + /i). The total superfluid density is ps = 7Pi + (1 — 7)p2 
with 7 = 1/(1 -|- fi). We also use the Debye temperature of 240 K [701 to calculate T^., 
which allows fixing one of the in-band pairing potentials. An. This leaves three free fit 
parameters: the second in-band potential, A22, inter-band coupling, A12, and the mass 
ratio, /i. Indeed, we found that Ps(T) can be well described in the entire temperature 
range by this clean-limit weak-coupling BCS model J5T]. Figure [2] shows the fit of 
the experimental superfluid density to the 7— model. The inserts show temperature - 
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dependent superconducting gaps obtained as a solution of the self-consistency equation, 



Eq. (34) and the lower inset show the gap ration as a function of temperature. Evidently, 



the smaller gap does not exhibit a BCS temperature dependence emphasizing the failure 
of the commonly used a— model. 

Lastly we note that in order to have two distinct gaps as observed in many 
experiments 0] one has to have significant in-band coupling constants. An and An. 
To support this conclusion we calculate the gap ratio 7a = Ai(0)/A2(0) for fixed A12 
and varying An and A22- Red lines in each graph show 7a = 2. For all cases, one needs 
significant A12 and varying An to reach this gap ratio. Therefore, the original simplified 
s± model with two identical Fermi surfaces and only interband coupling, A12 7^ and 
An = A22 = does not describe the experimentally found two distinct gaps in Fe-based 
superconductors ^71j . 




Figure 3. Gap ratio, 7a = Ai(0)/A2(0) for indicated fi (to check the effect of 
different partial densities of states) and A12. The magnitude of 7a is given by the 
intensity shown on the right. Red hnes in each graph show 7a = 2. 



4. Effects of Scattering 



The scattering by impurities strongly affects the London penetration depth even in the 
simplest case of non-magnetic impurities in materials with isotropic gap parameter and 
for scattering processes which can be characterized by the scalar (isotropic) scattering 
rate. The A anisotropy and magnetic impurities complicate the analysis of A(T), among 
other reasons due to the T^ suppression in these cases. For extended treatment the reader 
is referred to [S51 SD] and here we only summarize properties related to the discussion of 
our results. Introducing the non-magnetic scattering rate, 1/r, and magnetic scattering 
rate (spin - flip), l/r^, London penetration depth is expressed as: 



(A^ 



lQTT^e^TN{0){viVk) 



'ik 



E 



Jo 



A//o + h/2T- 



(41) 
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here 1/r^ = 1/r ± l/rm- Clearly, the evaluation of A(T) in the presence of magnetic 
impurities is quite involved. There is one limit, however, for which we have a simple 
analytic answer: 

Gapless limit: This is the case when r^ is close to 2/i/Ao(0), the critical value for which 
Tc = 0, i.e. TmA ^ 1. The resulting expression for the order parameter is remarkably 
simple 



A^ = 2n^{T^ -T^) (42) 

The result for superfluid density is valid in the entire temperature domain, T < T^. 



^2x-i 8TTe'^N{0){viVk) (. 2p„ p- , ^ 



where p = h/{27iTcT), p^ = h/{27iTcTjn) and p^ = p ± pm- For a short transport 
mean-free path p ^ p^ "we have Abrikosov-Gor'kov's result: 

C PPm 

The idea of strong pair-breaking and, perhaps, gapless superconductivity finds 
experimental evidence in Fe-based superconductors in form of scaling relations for the 
specific heat jump and a pre-factor of the quadratic temperature variation of A(T) 
[551110]. 

To summarize, in the case of superconductor with line nodes, impurity scattering 
will change linear temperature dependence of the penetration depth at T ^ Tc to 
become quadratic [73] . Disorder will also lift the c-axis line nodes in case of extended 
s— wave and induce a change from effective T^ to exponentially activated behavior for 
the in-plane penetration depth [71] . However, if we start in the clean limit of a fully - 
gapped superconductor and introduce pair-breaking scattering (either due to magnetic 
impurities or due to unconventional gap structure, such as s±) [71], penetration depth 
will also exhibit a power-law behavior approaching T^ variation in the gapless limit. 
Therefore, with increasing impurity scattering, if A(T) ~ T", we expect the exponent n 
to change from 1 to 2 in case of nodal superconductor (or even from 1 to ~ exp in case 
of non symmetry - imposed nodes) and from ~ exp to 2 in case of fully-gapped material 



with pair-breaking. This is schematically illustrated by shaded areas in Fig. 10 (where 
we used large values of the exponent n to designate ~ exp behavior). 

5. Experimental results 

Ba(Fei_2;Ta;)2As2 is one of the most studied systems among all Fe-based 
superconductors and we have collected extensive data that illustrate general features 
often common to many other members of the diverse pnictide family. Here we focus on 
the electron doped Ba(Fei_2;Ta,)2As2 with T = Co and Ni. One reason why these series 
were chosen is because large, high quality single crystals are available [TT]. All samples 
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were grown from the self-flux and were extensively characterized by transport, structural, 
thermal and magneto-optical analysis. They all exhibited uniform superconductivity 
at least at the 1 fim scale and dozens of samples were screened before entering into 
the resulting discussion [751 E!- To demonstrate sample quahty we show magneto- 
optical images in Fig. |4] and Fig. [5j Details of this visualization technique are described 
elsewhere [76]. In the images intensity is proportional to the local magnetic induction. 
All samples show excellent Meissner screening [70]. Fig. Ill shows penetration of the 
magnetic field into the optimally doped sample at 20 K. A distinct "Bean oblique wedge" 
shape [77] of the penetrating flux with the current turn angle of 45° (implying isotropic 
in-plane current density) is observed. 
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Figure 4. Magnetic flux penetration at 20 K into a crystal with x=0.07. The last 
frame shows a schematic overlay of the expected "Bean oblique wedge" shape with 
isotropic in-plane current density. 



Furthermore, to look for possible mesoscopic faults and inhomogeneities, we show 
the trapped magnetic flux obtained after cooling in a 1.5 kOe magnetic field and 
turning field off. The vortex distribution is quite homogeneous indicating robust uniform 
superconductivity for various doping levels. This is shown in Fig. |4] for four different 
doping levels. 

The parent compound, BaFe2As2 is a poor metal [TS] having a high temperature 
tetragonal phase with no long range magnetic order and undergoes structural and 
magnetic transitions around 140 K into a low temperature orthorhombic phase with 
long range antiferromagnetic spin density wave order [TT] . Transition metal doping onto 
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Figure 5. Magnetic flux trapped in samples after cooling in a 1.5 kOe magnetic field 
to 5 K and turning field off. Doping levels are indicated in right top corners. 



the iron site serves to suppress the structural and magnetic transition temperatures and 
superconductivity emerges after magnetism has sufficiently been weakened. Doping into 
barium site with potassium results in hole doped superconductivity. Properties of this 
hole - doped system, at least as far as penetration depth is concerned, are quite similar to 
the electron - doped FeT122 [32], [79]. On the other hand, properties materials obtained 
by isovalent doping of phosphorus into the arsenic site seems to induce superconductivity 
without introducing significant scattering and seems to result in a superconducting gap 
with line nodes [80] . 

In the charge - doped systems, angle-resolved spectroscopy (ARPES) [ST | I52 | [SH ] IM] 
and thermal conductivity [IHIHS] consistently show fully gapped Fermi surfaces (however 
with gap anisotropy increasing upon departure from optimal doping in case of thermal 
conductivity [IS].) Furthermore, the c— axis behavior is quite different and it is 
possible that a nodal state develops upon doping beyond optimal level [121 116]. The 
in-plane penetration depth consistently shows non-exponential power-law behavior 
[35l[8ail0l[86l|871|32l[88l|26l[27l[79l|89l[90ll^ which will be discussed in detail 

below. It seems that such behavior can be explained by the pair-breaking scattering 
[92| l93| [Ml [95], |96|. This is supported experimentally by the observed variation of the 
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low-temperature A(T) within nominally the same system (and even in pieces of the same 
sample) [32] as well as deliberately introduced defects [88]. In this review we also discuss 
the case of a substantial variation of A(T) between various samples, probably due to 
edge effect. 

In the following analysis we use two ways to represent the power law behavior: 
A(T) = A(0) + AiT/Tc)"^ at low temperatures (below 0.3 Tc) with A being the only free 
parameter, because at a gross level, all samples follow the A(T) ~ T^ behavior rather 
well and A(T) = A(0) + CT", leaving the exponent n as free parameter to analyze its 
evolution with doping or artificially introduced defects. In the case of vertical line nodes, 
we expect a variation from n = Ito n = 2 upon increase of pair-breaking scattering [75] , 
but in the case of a fully gapped s± state we expect an opposite trend to approach n = 2 
in the dirty limit from clean - limit exponential behavior [221 ESI EH ES]- If, however, 
nodes are formed predominantly along the c-axis in the extended s— wave scenario, the 
effect of scattering on the in-plane penetration depth would be opposite - starting from 
roughly n = 2 in the clean limit and approaching exponential in the dirty limit |96j . 

5.1. In-plane London penetration depth 

Figure |6] shows normalized differential TDR magnetic susceptibility of several single 
crystals of Ba(Fei_2;Coa;)2As2 across the superconducting "dome". All but one samples 
were grown at the same conditions and with similar starting chemicals. All, but one had 
thicknesses in the range of 100 - 400 nm. One of the samples {x = 0.074, denoted batch 
7^2) was cleaved for the irradiation experiments and had thickness of 20 nm. We use it 
for comparison with the "thick" batch #1 and, also, to study the effects of deliberately 
induced defects. It turns out that the edges of the thicker samples are not quite smooth 
and, when imaged in SEM, look hke a used book (see Fig. [ll|(a)). Since calibration of 
the TDR technique relies on the volume penetrated by the magnetic field, the thinner 
samples should be closer to the ideahzation of the sample geometry (top and bottom 
surfaces are always very flat and mirror-like), thus producing a more reliable calibration. 
On the other hand, this would only lead to a change of the amplitude (due to geometric 
mis-calibration) of the penetration depth variation (i.e., pre-factor A) and would not 
change its functional temperature dependence (i.e., the exponent n). Thinner samples, 
on the other hand, have better chance to be more chemically uniform, thus have reduced 
scattering. We observed these effects comparing samples from different batches. 

Low-temperature variation of London penetration depth is presented in Fig. [7] as 
a function of (T/T^)'' obtained by fitting the data to AA(T) = A(0) + C{T/T^Y. Each 
curve reveals a robust power law behavior with the exponent n shown in the inset in 
Fig. [7j The fitted exponent n varies from n = 2 ± 0.1 for underdoped samples to 
n = 2.5 ±0.1 for the overdoped samples within the batch 7^1 and reaches n = 2.83 
in batch 7^2. If the superconducting density itself follows a power law with a given 
n, then C ~ fs{c/ujp)S, where fs is the superconducting fraction at zero temperature, 
c is the speed of light and 5* is defined by the fraction of the Fermi surface that is 
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Figure 6. 47rx (T) in single crystals of Ba(Fei_a;Coa;)2As2 for different x. Sample 
from batch #1 is shown by a solid line. 



gapless (which may reflect a multigap character of the superconductivity, possible nodal 
structure, unitary impurity scattering strength, etc) and cOp is the plasma frequency. 
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Figure 7. Low-temperature behavior of AA(r) vs. {T/Tc}"' for all studied 
concentrations. Inset shows exponent n as function of concentration. Sample from 
batch =f/=2 is shown for comparison. 



Clearly the sample of batch ^2 shows behavior much closer to exponential compared 
to batch #1. As discussed above, this could be due to the variation of scattering between 
the batches. We analyze low-temperature A(T) in Fig. |8} We attempted to fit the data 
with three functions: the power-law with free pre-factor C and exponent n, the standard 
single-gap BCS behavior, Eq. (27), with a fixed value of A(0) = 200 nm and A(0) as a free 



parameter and to a BCS - like function where both A(0) and A(0) were free parameters. 
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The resulting values are shown in Fig. |8| The power-law fit yields quite high exponent 
n ^ 2.83 and the best fit quality. The BCS - like fit yields a reasonable fit quality, but 
produces impossible values of both A(0) ~ 48 nm and A(0) ~ 0.78Tc (the latter cannot 
be less than a weak coupling BCS value of 1.76. Finally, the fixed A(0) BCS fit does 
not agree with the data and also produces unreasonable A(0) ~ 1.25Tc. One strong 
conclusion follows from this exercise - we are dealing with a multi-gap superconductor. 
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Figure 8. Low-temperature behavior of AA(r) for tiie sample batch. #2. Solid lines 
are the fits to three functions described in the text. Insert shows full temperature 
variation indicating very sharp superconducting transition. 



In order to understand the validity of the empirical power-law behavior. Fig. |9] 
shows low-temperature behavior of AA(T) vs. (T/Tc)'^'^^ for the sample from batch 
7^2. Arrows show actual reduced temperature. Inset zooms at below of the commonly 
accepted "low-temperature limit" of ~ T^/S. Clearly, power-law behavior is robust and 
persists down to the lowest temperature of the experiment oi T ^ 0.02Tc. 

We now summarize the observed power-law behavior of the in-plane penetration 
depth for the electron - doped 122 family of superconductors. Figure [lO] shows the 



experimental low-temperature limit power-law exponent for different dopants on the 
Fe site and at different doping regimes. The shaded areas show the expectations 
for the pair-breaking scattering effects in s— and d— wave scenario. It seems that 
statistically rf— wave pairing (more generally - vertical line nodes) cannot explain the 
in-plane variation of the penetration depth. However, if the nodes appear somewhere 
predominantly along the c-axis, they may induce an apparent power law behavior of 
Xab{T) with the effective n f» 2 in the clean limit and going towards exponential behavior 
with the increase of the scattering rate [96] . 
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Figure 9. Low-temperature behavior of AA(r) vs. (T/Tc)'^-^^ for the sample from 
batch #2. Marks show actual reduced temperature. Insets show the behavior below 
commonly accepted "low-temperature limit" of sa T^/S. 
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Figure 10. Power-law exponent of the low-temperature variation of in-plane XatiT) is 
several electron-doped Bal22 superconductors at various doping levels. Shaded areas 
show the influence of pair-breaking scattering with n = 2 being the limiting value 
of n approaching from either the s— wave side (nodeless exponential) or d— wave side 
(vertical line nodes). In case of extended s— wave with nodes predominantly along the 
c-axis, effective n increasing from 2 in the clean limit towards the exponential behavior 
when the nodes are lifted by scattering. 



5.2. Absolute value of the penetration depth 

To further investigate the effects of doping and the difference between the batches, we 
use the method described in section 2.1. 2| which involves measuring the sample, coating 
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it with a uniform layer of Al and re-measuring [501 EH] . The Al film was deposited onto 
each sample while it was suspended from a rotating stage by a fine wire in an argon 
atmosphere of a magnetron sputtering system. Film thickness was checked using a 



scanning electron microscope in two ways, both of which are shown in Fig. 11 The first 
method involved breaking a coated sample after all measurements had been performed 
to expose its cross section. After this, it was mounted on an SEM sample holder using 
silver paste, shown in Fig. [ll](a). The images of the broken edge are shown for two 
different zoom levels in Fig. [TT|^b) and (c). The second method used a focused ion beam 
(FIB) to make a trench on the surface of a coated sample, with the trench depth being 
much greater than the Al coating thickness, shown in Fig. [TT|^d). The sample was then 
tilted and imaged by the SEM that is built into the FIB system, shown in Fig. fllle). 




Figure 11. Scanning electron microscope images of tlie Al coated samples, (a) Large 
scale view. The broken side is on top. (b) and (c) are zoomed in on the Al film on 
the edge of the broken side, (d) A trench produced by a focused-ion beam (FIB), (e) 
Close-up view of the FIB trench showing the Al film and its thickness. 

Example of the penetration depth measurements before and after coating are shown 
in Fig. [T2j Notice how small is the effect of Al coating when presented on a large scale 
of a full superconducting transition of the coated sample. However, TDR technique is 
well suited to resolve the variation due to aluminum layer [501 EH] . 



Obtained values of Aab(O) are summarized in the top panel of Fig. 13 for doping 



levels, X, across the superconducting region of the phase diagram, shown schematically 



in the bottom panel of Fig. 13 The size of the error bars for the Aafe(O) points was 
determined by considering the film thickness to be t = 100 ±10 nm and Xai{0) = 50 ±10 
nm. The scatter in the \ab{0) values shown in the upper panel of Fig. [13] has an 
approximately constant value of ± 0.075 yum for all values of x, which probably indicates 



that the source of the scatter is the same for all samples. For comparison. Fig. 13 also 
shows Xab{0) obtained from /xSR measurements (red stars) [22], the MFM technique 
(open stars) [251 EZ] and optical reflectivity (purple open triangles) [HZ] • Given statistical 
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Figure 12. Main frame : Full superconducting transition of an optimally doped 
FeCol22 crystal from batch #1 before and after Al coating. Inset : Zoomed in low- 
temperature region, T^m ^T < T^'' , before (green triangles) and after (brown circles) 
the Al coating on the same sample. (Notice how small is the effect on a full scale of 
the main frame). 



uncertainty these measurements are consistent with our results within the scatter. It 
may also be important to note that the Aa6(0) values from other experiments are all on 
the higher side of the scatter that exists within the TDR Aab(O) data set. As discussed 
above, our TDR techniques give a low bound of A(0), consistent with this observation. 
In order to provide a more quantitative explanation for the observed increase in 
Aab(O) as X decreases in the underdoped region, we have considered the case of s^ 
superconductivity coexisting with itinerant antiferromagnetism [SB]- For the case of 
particle hole symmetry (nested bands), the zero temperature value of the in-plane 
penetration depth in the region where the two phases coexist is given by 



A 



SC+SDW 
ab 



(0) 




(45) 



where A°^(0) is the value for a pure superconducting system with no magnetism 
present, and Aaf and Aq are the zero temperature values of the antiferromagnetic 
and superconducting gaps, respectively. Deviations from particle hole symmetry lead 



to a smaller increase in Af^''''^^^(0), making the result in Eq. (45) an upper estimate 

EHl. 



The three blue dashed lines shown in the top panel of Fig. 13 which were 



produced using Eq. (45), show the expected increase in Aab(O) in the region of coexisting 
phases below x ~ 0.047 by normalizing to three different values of Aa;,(0) in the pure 
superconducting state, with those being 120 nm, 180 nm and 270 nm to account for the 
quite large dispersion of the experimental values. This theory does not take into account 
changes in the pure superconducting state, so for x > 0.047 the dashed blue lines are 
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Figure 13. Top panel: The zero temperature London penetration deptti, Aab(O), as a 
function of tlic Co concentration, x. The three dashed blue lines are theoretical curves 



obtained using Eq. (45 1 for three different values of Xab{0) in the pure superconducting 
state. The solid gray line is a fit to the TDR data only of the form A+B/x". Also 
shown are values of AQf,(0) obtained by other experiments for comparison explained 
in the text. Bottom panel: Schematic phase diagram for FeT122 system showing the 
coexisting region [98, 99 . 



horizontal. These theoretical curves were produced using parameters that agree with 
the phase diagram in the bottom panel of Fig. 13 [HHl US]- While the exact functional 



form was not provided by any physical motivation and merely serves as a guide to the 



eye, the solid gray line in Fig. 13 is a fit of the TDR Aafc(O) data to a function of the 



form A+B/a;", which does indeed show a dramatic increase of Aab(O) in the coexistence 
region and also a relatively slight change in the pure superconducting phase. It should 
be noted that a dramatic increase in Aab(O) below x ~ 0.047 cannot be explained by the 
impurity scattering, which would only lead to relatively small corrections in A(T). 

With the experimental values of A(0) we can now analyze the superfluid density. 



In general A (T) is given by Eq. (23) and depends on the averaging over the 



particular Fermi surface. For example, in the simplest cylindrical case, Eq. (30) and 
A^^(O) = 47re^A^(0)f^/c^, where iV(0) is the density of states at the Fermi level and 
V is Fermi velocity. However, it is instructive to analyze the behavior from a two- 
fluid London theory point of view looking at the density of superconducting electrons. 
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Figure 14. Top panel: A~^ oc Ug vs. T for underdoped, optimally doped and 

overdoped concentrations for batch #1. A sample from batch #2 is shown for 

comparison. Bottom panel: normalized superfLuid density, ps = ns{T)/ns{0) = 
(A(0)/A(T))^ 



Us = (mc^/47re^)A~^ as function of temperature. The zero value, ns{0) will in general 
be less than total density of electrons due to pair-breaking scattering, so the magnitude 
of Hg is useful when comparing different samples. 



Figure [14] shows the data for underdoped, optimally doped and overdoped samples 
from batch #1 and also a sample from batch j^2 for comparison. For this sample #2 
A(0) = 200 nm was used. There is a clear, but expected, asymmetry with respect 
to doping. Underdoped samples show quite low density, because not all electrons are 
participating in forming the Cooper pairs and parts of the Fermi surface are gapped 
by the SDW as was discussed above. The overdoped sample, x = 0.01, despite having 
smaller Tc than the ones with x = 0.074, shows the highest n^. Obviously, the data 
scatter is significant. Therefore, the only reliable conclusion is that penetration depth 
increases dramatically upon entering the coexisting region. The overdoped side has 
to be studied more to acquire enough statistics. Furthermore, comparing two samples 
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with X = 0.074 from two different batches reveals an even more striking difference. Not 
only does the sample from batch #2 has larger n, but the temperature dependence 
of ps in the full temperature range is also quite different. The pronounces convex 
shape (positive curvature) of Ps(T) observed in all samples from batch #1 at the 
elevated temperatures becomes much less visible in sample 7^2. The bottom panel 



of [14] clearly demonstrates this difference, which is hard to understand based purely on 
the geometrical consideration (different thicknesses). It seems that thicker samples of 
the batch #1 had higher chance of being chemically inhomogeneous across the layers. 
On the other hand, the convex shape of Ps(T) at elevated temperatures is a sign of the 
two- gap superconductivity [56J, which depend sensitively on the interaction matrix, \^^, 



see section 3.2 and Eq. (34). This feature becomes more pronounced when the interband 
coupling, A12, becomes smaller compared to the in-band coupling potentials. An and 
A22- If our interpretation that the difference between x = 0.074 samples from batch #1 
and batch #2 is due to enhanced pair-breaking scattering in ^1, this would indicate 
that this scattering is primarily of inter-band character, so it disrupts the inter-band 
pairing. This gives indirect leverage to the s± scenario where inter-band coupling plays 
the major role. 
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Figure 15. Normalized superfluid density for sample ^2. Symbols show the data and 



the solid lines represent partial and total p{T) obtained from a fit to Eq.(40). Also 
shown are the clean (dashed grey lines) and dirty (dotted grey lines) single gap s— and 
d— wave cases. Inset shows superconducting gaps obtained self-consistently during the 



fitting, see Eq.(36l 



The normalized superfluid density for sample #2 is analyzed in Fig. 15 by using a 



two-band 7— model described is section |3.2[ Symbols show the data and the solid lines 
represent partial, pi(T) and P2(T), as well as total Ps(T) obtained from the fit using 



Eq.(40). The fit requires solving the self-consistent coupled "gap" equations, Eq. (36) 
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which are shown in the inset in Fig. [T5j Parameters of the fit are as follows: Ail = 0.80, 
A22 = 0.49, Ai2 = 0.061, 7 = 0.87. We used a Debye temperature of 250 K [4J to 
obtain the correct Tc via Eq. 37) that fixed An and gave A = 0.41. We also assumed 
equal partial densities of states on the two bands, so the value of 7 = 0.87 most likely 
comes from the difference in the k— dependent Fermi velocities, but may also reflect 
the fact that densities of states are not equal. Indeed, presented fitting parameters 
should not be taken too literally. The superfiuid density is calculated from the 



temperature - dependent gaps (inset in Fig. 15) and close temperature dependencies can 
be obtained with quite different fitting parameters. However, the gaps fully determine 
the experimental ps and this is the main result. We find that Ai(0) = 1.883Tc = 3.73 
meV and A2(0) = 0.754Tc = 1.49 meV, which are in good agreement with specific 
heat |100t [SH] and /iSR penetration depth measurements [H] done on the samples of 



similar composition. Also shown in Fig. 15 ) are the clean (dashed grey lines) and dirty 
(dotted grey lines) single gap s— and d— wave cases. (Note that while the gap does not 
depend on the non-magnetic impurities in isotropic s— wave case (Anderson theorem), 
the superfiuid density does.) Clearly, Ps(T) for sample #2 (and, of course for samples 
of batch #1) does not come even close to any of these single-gap scenarios. 
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Figure 16. In - plane, Xab{T), and out -of- plane, Ac(T) in a single crystal of FeCol22 
with X ~ 0.075. Inset - temperature - dependent anisotropy, 7^ — Xc/^ab- 



5.3. Anisotropy of London penetration depths 

Let us now discuss the electromagnetic anisotropy in the superconducting state, 
parameterized by the ratio 7a = Xc/^ab- To determine Ac we used method described 
in section 2.1.3| and the results are presented in Fig. 16 The problem is that we do 
not know the absolute value of Ac(0), so we could only obtain AAc(r) with the help of 
knowing AXab{T), which was measured on the same crystal independently. To find the 
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total Ac(T) we use the fact that close to Tc (in the region of validity of Ginzburg-Landau 
theory), we should have 



7A (TJ = ^7p (Te) (46) 

where anisotropy of normal state resistivity, 7p = Pc/pab is taken right above Tc. With 



7P ?a 4 ± 1 [lOlj . so that 'jxiTc) ~ 2. The results is show in the inset in Fig. 16 
Of course, there is some ambiguity in determining the exact anisotropy value, but the 
qualitative behavior does not change - the anisotropy increases upon cooling. With our 
estimate it reaches a modest value of 5 at low temperatures, which makes pnictides very 
different from high-T^ cuprates. This is opposite to a two-gap superconductor MgB2 
where 7a decreases upon cooling |102[ [65]. which may be due to different dimensionality 
of the Fermi sheets. In any case, temperature - dependent 7a can only arise in the case 
of a multi-gap superconductor. 

Next we examine the anisotropy of A(T) at different doping levels. This 
study was performed on FeNil22 samples and is reported in detail elsewhere |12]. 
Figure [l7|^a) summarizes the T{x) phase diagram showing structural/magnetic (T^m) 
and superconducting (Tc) transitions. The inset shows TDR measurements in a full 
temperature range for all concentrations used in this study. Figure [T7|^b) shows the 
low-temperature (T < 0.3Tc) behavior of Xab{T) for several Ni concentrations. The 
data plotted versus (T/Tc)^ are hnear for underdoped compositions and show a clear 
deviation towards a smaller power-law exponent (below temperatures marked by arrows 
in Fig. [iTJ^b)) for overdoped samples. While at moderate doping levels the results are 
fully consistent with our previous measurements in FeCo-122 [861 EZ], the behavior in the 
overdoped samples is clearly less quadratic. It should be noticed that in order to suppress 
Tc by the same amount, one needs a two times lower Ni concentration compared to Co. 
In previously FeCo-122 [86], the samples never reached highly overdoped compositions 
equivalent to x = 0.072 of Ni shown in Fig. flTFb). Therefore, Ni doping has the 
advantage of spanning the phase diagram with smaller concentrations of dopant ions, 
which may act as the scattering centers. The evolution of the exponent n with x is 



summarized in the lower inset in Fig. 17 



Now we apply a technique described in Sec. 2.1.3 to estimate Ac(T). Figure 18 shows 
the effective penetration depth, Xmix (see Eq. |8|, for overdoped, x=0.072 (main panel), 
and underdoped, x=0.033 (inset), samples before (A) and after (B) cutting in half along 
the longest side (/-side) as illustrated schematically at the top of the figure. Already 
in the raw data, it is apparent that the overdoped sample exhibits a much smaller 
exponent n compared to the \ab(T), while underdoped samples show a tendency to 
saturate below 0.13Tc. Using Eq. [8] we can now extract the true temperature dependent 



AAc(T). The result is shown in Fig. 19 for two different overdoped samples of the same 
composition, x = 0.072 having Tc = 7.5 K and T^ = 6.5 K, and for an underdoped 
sample with x = 0.033 having Tc = 15 K. Since the thickness of the sample is smaller 
than its width, we estimate the resolution of this procedure for AAc to be about 10 nm, 
which is much lower than 0.2 nm for AA^fe. Nevertheless, the difference between the 
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Figure 17. (a) Temperature-doping phase diagram of Ba(Fei_a;Nij;)2As2. The inset 
shows full-temperature range superconducting transitions, (b) AAaf,(r) for different 
doping levels vs. {T/Tc)"^. The arrows mark the temperature below which n becomes 
less than 2. The curves are shifted vertically for clarity. The inset shows the power-law 
exponent n{x) obtained by fitting to XatiT) = a + hT^ from the base temperature up 
to r/Tc=0.3. 



samples is obvious. The overdoped samples show a clear linear temperature variation 
up to Tc/3, strongly suggesting nodes in the superconducting gap. The average slope is 
large, about dX^/dT ^ 300 nm/K indicating a significant amount of thermally excited 
quasiparticles. By contrast, in the underdoped sample the inter-plane penetration depth 
saturates indicating a fully gapped state. If fitted to the power-law the exponent in the 
underdoped sample 2 < ra < 3, depending on the fitting range. 

Nodes, if present somewhere on the Fermi surface, affect the temperature 
dependence of both components of AA(T). However, the major contribution still comes 
from the direction of the supercurrent flow, thus placing the nodes in the present case 
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Figure 18. Schematics of magnetic field penetration in the case of Hac \\ I for the 
whole sample [A] and after cutting in half along length / [B] . (Main panel) The mixed 
penetration depth AXmix (T) before [A] and after [B] cutting for the overdoped sample 
x=0.072, Tc = 7.5 K. Inset shows similar data for the underdoped x=0.033, Tc = 15 
K. 



at or close to the poles of the Fermi surface. The nodal topologies that are consistent 
with our experimental results are latitudinal circular line nodes located at the finite 
kz wave vector or a point (or extended area) polar node with a nonlinear (A(0) ~ 6^, 
p > 1) variation of the superconducting gap with the polar angle, 9. It is interesting to 
note a close similarity to the results of thermal conductivity measurements in overdoped 
FeCol22 that have reached the same conclusions - in - plane state is anisotropic, but 
nodeless |15], whereas out of plane response is nodal [36]. Still, we emphasize that the 
apparent power-law behavior of the in-plane penetration depth, Xab{T), in a heavily 
overdoped samples could be induced by the out-of-plane nodes |lU3t l55|. To summarize, 
it appears that not only is the gap not universal across different pnictide families [80] . 
it is not universal even within the same family over an extended doping range. Similar 
conclusion has been reached for the hole-doped BaK122 pnictides |104t 11051 1106] . 



5.4- Pair-breaking 

Although the natural variation in the scattering rates between samples provides a good 
hint towards importance of pair-breaking scattering, for more quantitative conclusions 
we need to introduce additional disorder. This is can be achieved with the help of 
heavy-ion irradiation. To examine the effect of irradiation, ~ 2 x 0.5 x 0.02 — 0.05 mm^ 
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Figure 19. A\c{T) for the underdoped, x = 0.033, T^ = 15 K (UNDER), and for two 
ovcrdoped, x = 0.072, T^. = 7.5 K and T^, = 6.5 K (OVER), samples. Dashed hnes are 
hnear fits. 



single crystals were selected and then cut into several pieces preserving the width and 
the thickness. We compare sets of samples, where the samples in each set are parts 
of the same original large crystal and had identical temperature-dependent penetration 
depth in unirradiated state. (These samples is what we call batch #2 in this review 
with unirradiated reference piece appearing in the discussion and figures of the previous 
sections). Irradiation with 1.4 GeV 208p]^56+ j^j^g ^^ performed at the Argonne Tandem 
Linear Accelerator System (ATLAS) with an ion flux of ~ 5 x 10^^ ions-s~^-m~^. 
The actual total dose was recorded in each run. Such irradiation usually produces 
columnar defects or the elongated pockets of disturbed superconductivity along the ions 
propagation direction. The density of defects, d, per unit area is usually expressed in 
terms of so-called "matching field", 5^ = ^o*^, which is obtained assuming one flux 
quanta, $o ~ 2.07 x 10^'^ G-cm^ per ion track. Here we studied samples with B^ = 0.5, 



LO and 2.0 T corresponding to d = 2.4 x 10^° cm"^, 4.8 x 10^° cm"^ and 9.7 x 10^° cm"^. 
The sample thickness was chosen in the range of ~ 20 — 50/xm to be smaller than the 
ion penetration depth, 60 — 70 /xm. The same samples were studied by magneto-optical 
imaging. The strong Meissner screening and large uniform enhancement of pinning have 
shown that the irradiation has produced uniformly distributed defects |107j . 

Indeed, to see the effect, we need to start with the best (largest exponent) sample 
we have. We irradiated the sample designated as batch #2 with n = 2.83, which was 



discussed in detail above (see Figures |8| [9] and 15). To analyze the power-law behavior 
and its variation with irradiation, we plot AA as a function of {t = T/Tc)"'° in Fig. 20 



where the hq values for FeCol22 and FeNil22 were chosen from the best power-law 



fits of the unirradiated samples (see Fig. 21). While the data for unirradiated samples 



appear as almost perfect straight lines showing robust power-law behavior, the curves for 
irradiated samples show downturns at low temperatures indicating smaller exponents. 
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Figure 20. Detailed comparison of the functional form of AA(r) for irradiated 
FeCol22 and FeNil22. In the main panels AA(r) is plotted vs. (i = T/TJ"" with 
the exponents n taken from the best fits of unirradiated samples: hq — 2.8 and 2.5 for 
FeCol22 and FeNil22, respectively. Apparently, irradiation causes low-temperature 
deviations, which are better seen in the derivatives, dAX{t)/dt"", plotted in the insets. 



This observation, emphasized by the plots of the derivatives dAX{t)/dt"'° in the inset of 



Fig. 20 , points to a significant change in the low-energy excitations with radiation. 



The variations of Tc and n upon irradiation are illustrated in Fig. [21} Dashed lines 
and circles show FeCol22, while solid lines and triangles show FeNil22. The upper 
panel shows the variation of Tc and the width of the transition. Since Bip is directly 
proportional to the area density of the ions, d, we can say that T^ decreases roughly 
linearly with d. The same trend is evident for the exponent n shown in the lower panel 



of Fig. 21 



The influence of impurities, assuming s± pairing, has been analyzed numerically in 
a T-matrix approximation [88J • Figure [2l|^a) shows calculated superfluid densities for 



different values of the scattering rate. Figure 21 b) shows corresponding densities of 
states. Finally, Fig. 21 c) shows the central result: the correlation between Tc and n. 



Note that these two quantities are obtained independently of each other. Assuming that 
the unirradiated samples have some disorder due to doping, and scaling T"" to lie on 
the theoretical curve, we find that the Tc{B^) of the irradiated samples also follows this 
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Figure 21. Top panel: The suppression of Tc with disorder relative to unirradiated 
T"". The vertical bars denote the width of the transition corresponding the 
diamagnetic signal change from 90% (onset) to 20% (end). Symbols are shown at 
the mean values between onset and end of the transition. Lower panel: exponent n vs. 
Ba,. 



curve. The assumption of similarity between doping and radiation-induced disorder, 
implied in this comparison, while not unreasonable, deserves further scrutiny. More 
recent discussion on the variation of T^ with disorder is found in Ref. ( |108j ). 
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Figure 22. (a) Superfluid density and (b) the density of states, computed for the 
s± state with sign-changing isotropic gaps and strong interband impurity scattering, 
between the Born and unitary limits. The dashed line in (a) is an example of a power- 
law fit p{T)/po = p{0)/po - a(r/Tco)" for < T < 0.4 T^o; (b) As the impurity 
concentration riirnp ^ T)) increases, the band of mid-gap states approaches the Fermi 
level and the exponent n is reduced, (c) Tc vs. power n, from the theoretical model 
(triangles) and experiment (squares and circles). 
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6. Conclusions 

It was not possible to include in this review many interesting results obtained for 
various members of the diverse family of iron-based superconductors. However, we may 
provide some general conclusions based on our work as well as on results by others. 

(i) The superconducting gap in optimally doped pnictides is isotropic and nodeless. 

(ii) Pair-breaking scattering changes the clean - limit low temperature asymptotics 
(exponential for nodeless and T— linear for line nodes) to a ~ T^ behavior. 
Therefore, additional measurements (such as deliberately introduced disorder) are 
needed to make conclusions about the order parameter symmetry. In anisotropic 
superconductors in general and in s±, in particular, even the non-magnetic 
impurities are pair-breaking. 

(iii) The materials can be described within a self-consistent two-band 7— model with 
two gaps with the ratio of magnitudes of about Ai(0)/A2(0) ^2 — 3. 

(iv) Upon doping, the power-law exponent, n for the in-plane penetration depth, Xab{T), 
decreases reaching values below 2 signaling of developing significant anisotropy, 
whereas out-of-plane Ac(T) shows a linear — T behavior signaling of line nodes with 
Fermi velocity predominantly in the c— direction. 

(v) There is a fairly large region of coexisting superconductivity and long-range 
magnetic order, albeit with suppressed superfiuid density. 

(vi) Overall, the observed behavior is consistent with the s± pairing, but realistic 3D 
calculations are required to achieve agreement with experiments. 
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